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Abstract 

We prove that the motion of a triaxial Riemann elhpsoid of homogeneous 
hquid without angular momentum does not possess an additional first inte- 
gral which is meromorphic in position, impulsions, and the elliptic functions 
which appear in the potential, and thus is not integrable. We prove moreover 
that this system is not integrable even on a fixed energy level hypersurface. 

Keywords: Morales-Ramis theory. Elliptic functions, Monodromy, 
Differential Galois theory, Riemann surfaces 

1. Introduction 

We want to study the integrability of the following Hamiltonian 

oo 

2 , ^2^2 \ , f « 



with r = a/^i + g| (a notation that will always be used in the following). This 
Hamiltonian appears through the Hamiltonian formulation of the motion of a 
self gravitating triaxial Riemann ellipsoid of homogeneous liquid, in which we 
have restricted our study to the zero angular momentum case. This equation 
can be found in [l[ . The motion of Riemann ellipsoids have been studied for a 
long time, as for example in j2|,0]- Even the integrability of the Hamiltonian 
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H has been studied by Ziglin in This corresponds in fact to the hmiting 
case when we take an energy level which tends to infinity. However, this does 
not prove "meromorphic non integrability" of in a reasonable sense (nor 
gives a reasonable notion of "meromorphic non integrability" for H, as H is 
not meromorphic in p, q itself). This is because the first integral should have 
a "limit" when energy tend to infinity, and that this limit is meromorphic. 
It would have a "limit" if the first integral be rational in pi,P2, but even in 
this case, the question of the meromorphic dependence in the positions qi, q2 
would stay open. As the case a = has been already treated in j3| , one will 
only consider the case a ^ 0, which comes down after variable change to 
a = 1. 

The Hamiltonian H has two degrees of freedom, and there is an elliptic 
function in the potential. Because of this, one should precise the notion of 
integrability we want to study before any non integrability proof. Meromor- 
phic functions is clearly not a large enough class because it does not contain 
the Hamiltonian itself. So, we will not use directly Morales Ramis theory, but 
begin with definitions of what could be a reasonable meaning of integrability 
of H. Remark that this Hamiltonian is in fact well defined on a Riemann 
surface because of the elliptic integral, which is multivalued for (gi, ^2) ^ C^. 
This is here the main difficulty because contrary to algebraic functions, this 
transcendental extension cannot be suppressed in some additional first in- 
tegral only by algebraic transformations. So, a little expansion of Morales 
Ramis theory will be needed to deal with such Hamiltonians. 

Definition 1. We pose 

Jill, 12)= / , =dz 

and we will note S( J) the singular set of J, which is 

= {(gi,g2) e C^ J is not C°° on (51,^2)} 

Definition 2. We pose 

K = {g (gi, g2, Pi,P2, r, J, dg^ J, dg^ J) , g meromorphic} 

the (differential) field of meromorphic functions in p, q, r, J. Moreover, we 
will say 
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• that J G -ft' is a meromorphic first integral of if J is constant along 
any orbit. 

• that / G K is a meromorphic first integral on the level H = E if I is 
constant along any orbit with energy H = E. 

• that / G if is an additional first integral of H, or only on H = E, 
if I is functionally independent with H, or non constant on H = E 
respectively. 

Remark 1. Remark that the field of meromorphic functions 7^ is a differ- 
ential field, and contains the Hamiltonian. Remark also that even if a first 
integral on a level H = E is defined a priori everywhere according to Defini- 
tion [21 in fact we only care about its restriction on H = E. 

The "tool" we will use to prove non integrability is the Morales Ramis 
theory in [5|,[6|. 

Theorem 1. (Morales Ramis [^J Theorem 4-1- page 85) Let H he a Hamil- 
tonian, r a particular (not a point) orbit and suppose H meromorphic on 
a neighbourhood of T . If H possess a complete system of first integrals in 
involution, functionally independent and meromorphic on a neigbourhood of 
T , then the identity component of the Galois group of variational equation is 
abelian at any order. 

If the Hamiltonian satisfy good properties (as being meromorphic along 
the curve F), then a meromorphic first integral of H will produce a mero- 
morphic initial form on F and then a constraint on the Galois group. This is 
a very general way to prove non integrability: If a Hamiltonian is integrable 
in some sense, then the Galois group of the variational equation near a par- 
ticular orbit should satisfy some particular property. In their article, they 
consider meromorphic Hamiltonians and first integrals, and here we want to 
add elliptic functions, and prove that Morales Ramis approach is still valid 
in this case. 

The main theorem of this article is the following 

Theorem 2. The Hamiltonian H does not possess a first integral in K func- 
tionally independent with the Hamiltonian H . On a fixed energy hypersurface 
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H — E, the Hamiltonian system restricted to this hypersurface does not pos- 
sess a non constant first integral in K. 

The plan of the proof is the following 

• First we study smoothness of possible first integrals in in particular 
singular points of elliptic functions involved. 

• We compute a family of particular solutions corresponding to the case 
where the ellipsoid of fluid is invariant by rotation. A closed form so- 
lution can be given in this case, and we analyze (complex) singularities 
of these solutions. 

• We compute the variational equation near these orbits and prove a 
small extension of Morales Ramis for proving the non existence of an 
additional first integral in K. 

• As the variational equation has transcendental functions in its coeffi- 
cients, Kovacic algorithm cannot be used directly. Instead, we prove 
using an analysis of the monodromy that some second order differential 
equation related to the variational equation should have a Galois group 
whose identity component is solvable (which is a weaker property than 
the classical virtual abehanity condition from Morales Ramis) . 

By chance, even if the Hamiltonian has elliptic functions, in the computation 
we will never see anything "worse" than that. The difficulty is that comput- 
ing Galois group with transcendental functions in the base field can be much 
more tricky than in the rational base field case. 

2. Variational equation 

2.1. Regularity of function J 

We want to consider first integrals that will involve the function J, so we 
will need to know where they arc mcromorphic in a neighbourhood of qi = 0, 
which is the plane we will be interested in. We get the following lemma. 

Lemma 3. We have 

E( J) n {gi = 0} C {(0, 0), (O, 2ef ) , . . . , (o, 26^^) } 
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Proof. First we remark that J is an elliptic integral, and so is smooth in its 
parameters for values of the parameters for which the corresponding lattice 
does not collapse. This lattice collapses when the associated elliptic curve 
becomes singular, and it corresponds to the case where the polynomial of 
degree 3 under the square root has a double root. Here we have 



2 I I -'- „2 



z + ^ \ [ z +rz + -go 



gi=0 ^^2 



-^izql + A)i2z±q2)' 



So we cannot conclude directly. The function J can be written using classical 
elliptic integrals 



EUipticK ( ^ /zxTCT 1 - EllipticF ' / — ^ — / 21 



The denominator on gi = is equal to 



±ql + 8 

and so is smooth outside the seven points we exclude in Lemma [3l Further- 
more, the functions EllipticK(z), EllipticF (w, 2;) are even in z, and the sec- 
ond is smooth except for w = ±1. So the part in EllipticK is always smooth 
(outside the exceptional points), and the function EllipticF is smooth outside 
maybe the points such 

g2 (r - gi) 

For gi = 0, this constraint comes down to gf = ±8, the same constraint as 
before. So 

S( J) n {gi = 0} C {(0, 0), (0, 2ef ) ,...,(0,26"^)} 

2.2. The plane (pi,gi) = (0,0) 

To study integrability, we will need to look at a particular (explicit) so- 
lution. The solution we will study corresponds in fact to the case where the 
ellipsoid of fluid is invariant by rotation. 
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Theorem 4. The plane {pi, qi) = (0, 0) is an invariant plane. On this plane, 
the Hamiltonian can be written 

„5„2 arccos ( 2^/2,0^^^"^ 



+ <l2 ^Jql - 8 

This function is well defined on a (transcendental) Riemann surface outside 
its singularities which correspond to 



2i7v 4z7r 2i7r 4z7r 



q2 = 0, —1, — e a ,—6 3 , 2, 2e 3 , 2e 3 

On a each level R = E, there is only one orbit except for the energy values 
such that 

o AT 3/2 6V2V8^\ ^ , „ 12x 
2V2 - x'^' cosh = and E 



y x^ + IQ I x'^ + 16 

for which there is at least one fixed point and at most finitely many on the 
Riemann surface associated to R. 

Proof. The function R can be found by direct computation. There is a 
choice of valuation for the square root r, and when we take gi = this 
leads to terms of the type \/q^. We can choose in principle any of the two 
valuations (but we need to choose always the same in all the computations), 
and here we will choose -\/^ = q2- The fact that the elliptic integral can be 
now written with an arccosinus is not surprising because we already found in 
Lemma [3] that the corresponding elliptic curve becomes singular for qi = 0. 
Now this function R defines a one degree of freedom Hamiltonian, and so the 
orbit can be completely studied analyzing the levels of R. A critical point 
always corresponds to p2 = 0. We derive R in q2 for p2 = 

iV2{ql + 16) arccos {2^/2qo^'^\ ^ 

\ J ^ (1) 



2(8-g3)3/2 8-gf 
Then, substituting the arccos term using the equation R = E, we find 

gi + 16 



6 



This gives the equation of Theorem |H Now we check that such a critical 
point cannot be critical for 2 sheaves of the Riemann surface associated to 
R. If it was the case, then the equation would vanish for at least two 
valuations of the function arccos. Then we would have + 16 = 0. We 
check that these points are not critical points. So there are at most 3 critical 
points. 



We have in particular the level E = 1 with q2 = 2, which is an already 
known real fixed point, the others being complex. Typically, these critical 
levels lead to confluences in the variational equation and so have fewer con- 
ditions for integrability. However, we will see that in our case, this will have 
no impact on the non integrability proof. 

2.3. Normal variational equation 

Theorem 5. The normal variational equation near a (non constant) orbit 
in the plane (pi, qi) = with R = E is given by 

((4t^^ + 192t^ - 60t^ + 256t^)s - SOt'^E + 96t^E + 2t^^E + 12UE)X+ 

((75t^ - 640t - 7t^° - 72t^)s 
-384E + 42t^ - ge^t^ + 48t + 42t^E - ?,Et^ - 6t^)X+ 



(2) 



with 



((-64r + 2t'')s - 64r - 4r)X = 



. arccos 



^2^8^ 



Proof. We begin by direct computations of the Hessian matrix of H for 
= (0,0). This works in all cases except for dq^q^H.dq^q^H . We have 
in particular a variational equation of the form 



Y 



1 
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So the variational equation is already decoupled, and the normal part corre- 
sponds to the first 2x2 block. We have noted here Fx = dq^q^H. Let us look 
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now closer to this function. The kinetic part after derivation disappear for 
(Pi) <?!) = (0, 0). So the only thing to compute is 



' ^zq^ + 4{q2 + 2zf 



This integral can be computed without elliptic functions, only using the 
arccos function. After simplification, we find 

iq2\/2{qi — 32) arccos (^2\/2q2 j iq _|_ ^3 



4(8 - (8 - gi)2g2 

We now transform the normal variational equation in a second order differ- 
ential equation. This gives 

where {p2(t) , q2{t)) = is an orbit of H in the plane = 0. 

We now make the variable change — > t. For this, we use the fact that 
R is constant over an orbit and we get the relation 

. 2 ^72(03 + 1) arccos (2^20-3/2) ^3 ^ ^ 

(t> = , h E — 

03^/8303 04 

and also a relation for by derivating it. The variable change 0(t) — > t 
involves only and 0^. This eventually produces equation ([2]). 



3. Non-integrability 

Theorem 6. If H has an additional meromorphic first integral, then the 
Galois group over the base (differential) field K generated by meromorphic 
functions in 

\rt, V8 - arccos (l.^Plt-^l'^^ for t{f - + 1) ^ 

of equation ([2]) is virtually abelian for any E E C If H has an additional 
meromorphic first integral on some fixed energy level H = E, then the Galois 
group over K of equation ([2]) is virtually abelian for the corresponding value 
ofE. 
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Proof. Such a theorem is in fact a small extension of the Morales Ramis 
Theorem [TJ Suppose there is an additional meromorphic first integral / (in 
the sense of Definition [2]). We consider the curve given by 

r^: Pi = 0, qi = 0, R{p2,q2) = E, (0, gs) ^ 

excluding the possible critical points of R. The points (0,52) G S(-f^) cor- 
respond to singular points of the Hamiltonian, and thus we exclude them 
(there are finitely many because of Lemma [3]) . Remark that as proven in 
Theorem m the curve Te always exists (it is never only a point), there is a 
unique choice for Ve in R{p2,q2) = E and there is a discrete set of points 
removed from Te (recall that F^; is in the complex a Riemann surface, so 
points of S(if) can correspond to infinitely many points on F^;, one copy 
on each sheave). Let us prove that a first integral I & K is meromorphic 
along F^;. The first integral / can be written as a meromorphic function of 
p, q, r, J, dqJ. On Te, r has a singularity at g = 0, but is is also a singularity 
of J, and so is not in F^;. The singularities of J are not in F^; by definition, 
and dqJ are also elliptic functions whose singularities are at the same points 
as J. So we can apply Theorem [H We now need to check that the base 
field given in Theorem |6] is the good one, meaning the field of meromorphic 
functions on F^;. This corresponds to functions in P2,q2, meromorphic on 
a neighbourhood of the curve R{p2, ^2) = E outside the singular points, re- 
stricted to the curve R{p2i Q'2) = E. Using the same parametrization as in 
Theorem |5l we express p2 in function of q2 = t and this gives 



; i^/2{t^ + 1) arccos (2^21-^^) + 1 

P2 = \ , + E — - — 

After simplification, we see that the corresponding field is an algebraic ex- 
tension of degree 2 of and so the identity component of the Galois group 
of equation ([2]) will not change if we take K instead of the field as base field. 

The case of fixed energy is proven by reducing the dynamical system 
to the hypersurface H = E. We loose (at least a priori) the Hamiltonian 
structure and the notion of involution, but it was not needed because we are 
in two degrees of freedom (a presentation of the problem arising with more 
degree of freedom is done in [7|). The normal variational equation of this 3 
dimensional system is still equation (|2]), and the previous computations are 
still valid (we already proved that there are no "exceptional" energy levels, 
meaning that the curve F^; always exists and is never too singular). 
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Remark 2. Remark that the field K corresponds also to the field of mero- 
morphic functions on the Riemann surface 

S = {{x,y) eC^, x^cos^y = 8, - 8)(a;^ + 1) ^ 0} 

(replacing t by a; in the generators of K produce meromorphic functions in 
x,y, and reciprocally). We also have that the coefficients of equation (|2]) are 
meromorphic functions on S. In the following, we will keep this notation for 
the Riemann surface S. 

Theorem 7. Suppose that equation (j2]) has a virtually abelian Galois group 
over the base field K (for some fixed energy E). We have the following 

i The monodromy group Gi = {-R^, 7 closed curve on S} (R^ is the re- 
solvant matrix of equation is such that 

WgeV^'\G^), g'' = id 

where T>^'^\G) is the second derivative of the group G. 

ii The monodromy group G2 = {-R7, 7 closed curve on S} where is the 
resolvant matrix of the equation 

(4t" + 192t^ - 60t® + 2hQt^)X + (75t^ - 640t - 7t^° - 72t^)X+ 

(-64t^ + 2t^^)X = ^ 

is such that 

WgeV^'\G2), g'' = id 

iii The identity component of Galois group of ([3]) over the base field C(t) is 
solvable. 

Proof. Let us prove {%). Because the base field for Galois group compu- 
tations is not C(t), we cannot use Kovacic algorithm (at least not directly). 
But the classification of possible Galois group still holds. So the possible 
Galois groups are triangular groups, Dqo, finite groups D„, A4, 5*4, A^. Let us 
take one of these possible groups and derive it two times (the derivation of a 
group correspond to the group generated by its commutators). For triangu- 
lar groups, Daa, Dn, the second derivative produce the identity group. We 
also have 

V^^\Ai) = id V^^\Si) = TL\ V^'^\A^) = A5 
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So, in all cases, the elements of the second derivative of the Galois group have 
always an order dividing 60. We conclude using the fact that the monodromy 
group Gi is always a subgroup of the Galois group. 
Let us prove (u). We consider the application 



a 



: K — y K cr ^arccos = arccos 



The application a corresponds in fact only a translation of sheaves (x, y) — > 
{x,y + 27r) on the Riemann surface S. We now apply cr'^ on equation ([2]). 
This produces the equation 

( (At'' + 192t5 - 60t« + 256t2) ( s + 

-SOfE + m^E + 2t'^E + 12UE) X+ 
((75<'-640<-7<.«-72<-)(. + -^) (4) 

+ 42t^ - 96^t^ + 48t + A2t^E - 3Et^ - 6f) X+ 

( (-64t* + 2t'') (s + \ _ 54^5 _ \ ^ ^ Q 

^ V V2V8^J J 

So equation (jl]) will also have a first integral with coefficients in K. We 
consider the polynomial P = Az'' + 1922;^ — 60^;^ + 256^^, u > 0,e > two 
real numbers and the compact 

= {{D{0,iy) \ (u,eP-i(o)/^(^,e))) x D{0,iy))nS 

where D{z,e) is a disk with center z of radius e. For k > large enough, 
the singularities of equation ([2]) are not in because for large fc, the 
singularities of equation (j2]) inside Z)(0, z/)^ fl 5 are converging to the roots 
of P and — 8 = (and singular points of S are also roots of P). Now let 
consider four paths 71, . . . , 74 on 5 outside the roots of P and the commutator 

[[-R71 ' -R72] > [-R73 ) -^74]]^° (5) 

where i?-y is the resolvant matrix of dl]) along 7. There exist z/, e such that 
7i,...,74 C Cy^^. We divide (jl]) by its dominant term and we take the 
limit k — y 00. The equation is then converging to the equation ([3]). The 
resolvant matrix R is smooth when k — > 00 on the compact set Cu^e- So, for 
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k large enough, the resolvant matrices R are smooth along 7^, and so they 
are converging to monodromy matrices along curves 7^ of equation So is 
the commutator ([5]), which is equal to id by hypothesis. So 

^geV^'\G2), g'' = id 

Let us prove {Hi). We consider the monodromy group 

G3 = {R^, 7 closed curve on C} 

where Rj is the resolvant matrix of the equation ([3]). This monodromy group 
is a priori bigger than G2, because of multivaluated functions in K. We have 
that 

^ ( arccos ( 2^/21-""^^)) = , ^ o{C{t\, K) = 

dt\ V ^t5(8-t3) V ; 

Noting F the Picard Vessiot field of equation ([3]) and Kq the field of mero- 
morphic functions on {x G C, x{x^ — 8)(x^ + 1) 7^ 0}, we have a tower of 
fields Kq d K d F and monodromy groups 

a{K^,k) = D^ a{K,F)=G2 a{Ko,F) = Gs 

So, to "kill" this D^o that could appear in G3, we just need to derivate one 
more time. We get then 

ygeV^'HG,), g''' = id 

where we take the power 120 instead 60 because V{Doo) = Z2, and so the 
additional derivation does not kill it completely. Analysis of singularities of 
equation ([3]) shows that all its singularities are regular, and then that the 
equation is Fuchsian. This implies that the Galois group of equation over 
C{t) is the Zariski closure of the monodromy group G3. Suppose now that 
the Galois group of equation ([3]) over C(t) has not its identity component 
solvable. Following Kovacic's algorithm, the only possible groups contain 
(C). So one would need that G3 is Zariski dense in SL2{C). We have 

that 

V^'\SL2iC)) = SL2iC) ygeSL2iC), 3^ G ^^^(C), -g'""^ = g 

The application commutator and elevating to power 120 are continuous (in 
the Zariski topology, rational fractions are continuous), so this would imply 
that id is dense in SL2{C). Impossible. 
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Let us conclude by proving Theorem [21 

Proof. If H has an additional first integral in K (respectively on some 
particular energy level H = E), then using Theorem [6l we have that equa- 
tion ([2]) has a virtually abelian Galois group over the base field K, and then 
with Theorem [7] that equation (|3]) has a Galois group whose identity com- 
ponent is solvable. We just apply Kovacic's algorithm on equation ([3]), and 
find that its Galois group is SL2{C), and so its identity component is not 
solvable. This gives Theorem |2l 

4. Conclusion 

It is not so rare that a transcendental function appear in the variational 
equation we need to study to prove non integrability, but in general it is 
possible to avoid it by just taking a particular orbit for which this case does 
not occur, as done for example in |8i]. But we see that in fact it probably 
produces even stronger integrability conditions. Such a study is not so much 
more difficult because through a limiting process, such an equation comes 
down to an equation with coefficients in C(t) (the difference solvable abelian 
does not change anything in practice). And moreover, it is only a necessary 
criterion, if it was met, we could produce other conditions by making an 
asymptotic expansion in the "multivaluation parameter" (which corresponds 
here to apply the translation of sheave a). In the case where these transcen- 
dental extensions are not avoidable, this approach completely make sense 
because we do not have an equivalent of Kovacic algorithm for these. In the 
same problem of Riemann ellipsoid motion, zero angular momentum is only 
one case. A (probably) complete list of integrable cases is given in and 
it could maybe be possible to prove the non integrability of the other cases 
using this approach. 
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